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Abstract 

The ADM approach to canonical general relativity combined with Dirac's 
method of quantizing constrained systems leads to the Wheeler-DeWitt equa- 
tion. A number of mathematical as well as physical difficulties that arise 
in connection with this equation may be circumvented if one employs a co- 
variant Hamiltonian method in conjunction with a recently developed, math- 
ematically rigorous technique to quantize constrained systems using Rieffel 
induction. The classical constraints are cleanly separated into four compo- 
nents of a covariant momentum map coming from the diffeomorphism group 
of spacetime, each of which is linear in the canonical momenta, plus a single 
finite-dimensional quadratic constraint that arises in any theory, parametrized 
or not. 

The new quantization method is carried through in a minisuperspace ex- 
ample, and is found to produce a "wavefunction of the universe" . This differs 
from the proposals of both Vilenkin and Hartle-Hawking for a closed FRW 
universe, but happens to coincide with the latter in the open case. 



*E.P.S.R.C. Advanced Research Fellow 



1 Introduction 



Dirac's theory of constrained systems in classical mechanics (cf. Gotay et.al. (1978) 
for a modern geometric formulation) consists of two steps. Firstly, the constraints 
$ a = are imposed on the phase space S of the unconstrained system, singling out 
the constraint hypersurface C C S as their solution space. Secondly, one forms the 
quotient S* = C/JF of C by the foliation JF defined by the null directions of the 
induced symplectic form i*u on S (here ui is the symplectic form on S and i is the 
injection of C into S). This second step (which is absent if all constraints are second 
class) identifies physically equivalent points on C . If S — T*Q is a cotangent bundle 
and the constraints are components of a momentum map $ (cf. Marsden and Ratiu 
1994) derived from an action of a group G on Q, then all <£> a are necessarily linear in 
the momenta. The so-called Marsden- Weinstein quotient $ _1 (0)/G then coincides 
with the reduced phase space. 

In canonical (ADM) classical gravity (cf. Fischer and Marsden 1979, Isham 1993, 
Kuchar 1992) the configuration space Q is taken to be the space of Riemannian 3- 
metrics (subject to certain regularity conditions) on a (Cauchy) hypersurface S (here 
assumed to be compact) in space-time X. The $ a are the (super) Hamiltonian- and 
momentum constraints (one for each point x of S; a = 0, 1,2,3), the first of 
which are quadratic in the canonical momenta. Also, the Poisson bracket of two 
Hamiltonian constraints is proportional to the inverse of the three-metric, which 
makes it impossible to formulate these constraints as the pull-back of a Poisson 
morphism $ from S into any Poisson manifold P. (The super-momentum con- 
straints, which are linear in the canonical momenta, are of such a nature, with P 
being the dual of the Lie algebra of Diff(S) and $ the equivariant momentum map 
corresponding to the natural action of this group on S.) In particular, the connec- 
tion between the Hamiltonian constraints and Diff(X) is obscure. Time-evolution 
is generated by the constraints smeared with arbitrary functions iV(x, t); hence the 
dynamical evolution takes place along some of the null directions of and col- 
lapses to no evolution whatsoever on the physical phase space S°. With obvious 
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modifications, this discussion applies to the Ashtekar variables as well. 

In quantizing constrained systems, the two-step classcial procedure is replaced 
by a single step. Dirac (1964) singled out the first step: if Tt is the Hilbert space 
of states of the unconstrained system (that is, the quantization of S) and $ a are 
self-adjoint operators on H quantizing the classical constraints, then the quantum 
analogue Hd of the physical phase space S° is defined as Hd = H\^ a ip = Va}. 
This space then inherits the inner product from H, and is a Hilbert space in its own 
right. In other words, one imposes the constraints and that's it. This only makes 
sense if all constraints have as discrete eigenvalue, with common eigenspace. This 
condition is rarely satisfied in practice, and this has led to certain modifications of 
the Dirac proposal (Ashtekar and Tate 1994, Hajicek 1994, Ashtekar et.al. 1995, 
Marolf 1995b), in which one solves the constraints on a bigger space than H. One 
thereby loses the inner product, and the main problem is then to construct an inner 
product on the solution space from scratch. In some examples, such methods lead 
to the same result as our approach, but in the former one still tries to mimick the 
first step of the classical reduction process. Moreover, it is easy to devise examples 
where such rigged Hilbert space techniques will not lead to the desired answer, see 
below. 

Applied to canonical gravity (Ehlers and Friedrich 1994), the Dirac procedure 
leads to the Wheeler-DeWitt equation, in which the quantized Hamiltonian con- 
straints are imposed on the wavefunction ijj, which is a function(al) of the three- 
metric on E (or of the corresponding Ashtekar variables). Quite apart from the 
fact to what extent H is well-defined (for considerable progress in this direction see 
Ashtekar et.al. 1994, 1995), this approach meets formidable obstacles (Kuchar 1992, 
Isham 1993). 

The problems with the Wheeler-DeWitt equation can be traced back to (at least) 
two sources: i) the lack of covariance of the ADM approach, which is especially 
dangerous in connection with quantum field theory, where fields are not defined at 
sharp times; ii) the Dirac quantization method of constrained systems. It seems that 
one can do better on both accounts. Firstly, there exists a covariant Hamiltonian 
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formulation of classical field theory (Kijowksi 1973, Gotay 1991, Gotay et.al. 1993), 
and secondly the author (Landsman 1995) has recently formulated a new method of 
quantizing constrained systems, which avoids many problems in the Dirac approach. 

2 Constrained quantization revisited 

To start with the latter, the main idea is to implement a quantized version of the 
second step of the classical constraint procedure, rather than the first one. The con- 
straints then do not have to be implemented, avoiding the Wheeler-DeWitt equation 
altogether. This is done by modifying the inner product ( , ) on Tl (which is posi- 
tive definite) into a positive semi-definite sesquilinear form ( , )q. The latter has a 
nontrivial null space Afo = {ip G H\(ip,ip)o = 0}, in terms of which the state space 
H° of the constrained system is simply given by H° = H/Mq. Crucially, H° inherits 
the form ( , ) , which is now positive definite since its null vectors have been thrown 
away. Thus one has a bona fide inner product on TC°, which can be used to calculate 
physical amplitudes and probabilities. 

To illustrate the method, consider a compact group G, which acts on H through 
a unitary representation U . If the constraints correspond to the generators T a of (the 
Lie algebra of) G, the modified inner product turns out to be given by (Landsman 
1995) (■?/>, y?)o — fGdg{U(g)ip,<p). One recognizes that the right-hand side coincides 
with (P ip , P {p) , where P is the orthogonal projector onto the subspace 7i of 7i 
which transforms trivially under G. Hence Af = H^, so that H° = H , which 
coincides with Hd of Dirac. The reason one can get away with not explicitly imple- 
menting the constraints is that states not satisfying them are of zero norm in the 
modified inner product, and therefore quotient to the null vector in H . Also, the 
action of a constraint on a vector maps it into a null vector (in other words, a gauge 
transformation changes a vector only by a null vector). 

Had G not been compact, the modified form would be ill-defined on some of 
H, but our method is easily adapted to such instances: instead of all of H one 
finds a dense subspace L C H on which ( , ) is defined. The only change in the 
procedure is then that the quotient L/Af is not complete under the induced norm, 
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and has to be completed to form Tt . In this way, constraints with continuous 
spectrum can be handled without any difficulty (Landsman 1995). The general way 
of finding the modified inner product is based on the mathematical analogy between 
the momentum map and generalized Marsden-Weinstein reduction in symplectic 
geometry, and Hilbert C*-modules and Rieffel induction (Rieffel 1974) in operator 
algebra theory (cf. Landsman 1995 for details and references). 

For an interesting example which highlights the difference between our Rief- 
fel induction method and rigged Hilbert space techniques, take G a noncompact 
semisimple Lie group, and S = T*G with momentum map $ defined by the natural 
right-action of G on S. The classical reduced space $~ 1 (0)/G' is a point. Quanti- 
zation with our method proceeds by taking H = L 2 {G) as the quantization of S, 
which carries the right-regular representation U of G. Let ( , )o on L — C C (G) be 
defined as in the compact case. The physical Hilbert space comes out as H° = C, 
which is obviously the correct answer. However, the trivial representation of G is 
not weakly contained in the regular representation on L 2 {G), so using the standard 
rigged Hilbert space S(G) C L 2 (G) C S(G)' , where S(G) is the Harish-Chandra 
Schwartz space of G, would produce an empty physical state space. (Note that the 
word 'rigged' as used in Landsman (1995) to describe the modified inner product 
has nothing to do with the same word occurring in 'rigged Hilbert space'.) 

Let Vip G H° be the image of ip G L in the quotient L/Af , so that V : L — > H° 
satisfies (Vip,Vip) n o = (ip,(p) . In our examples, H° is of the form L 2 (R n ) (in 
momentum space), and V has the form (Vif>)(p) = (if), f p ) for a family of functions 
f p lying in some appropriate dual of L; the r.h.s. can be calculated as if it were the 
inner product in H, although f p is not in H, cf. Poerschke and Stolz 1993. 

Observables are those self-adjoint operators A on H which satisfy (Aif>, (p) = 
(if), A<p)o for all if), tp in L; this condition is the quantum version of the classical 
condition that {/, $ a } vanish on C. In either case, the point is that observables are 
well-defined on the reduced spaces Ti° or S* , respectively. The physical action tt° 
on H° is given by n°(A)Vip = VAif). 
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3 Covariant Hamiltonian method 



For our quantization method to apply, the classical reduced space has to have the 
form of a generalized Marsden-Weinstein quotient. This renders our method inap- 
plicable to ADM gravity. However, the covariant Hamiltonian approach advertised 
above creates a context in which our technique does apply. For the parametrized 
particle and minisuperspace examples we will have to deal with noncompact abelian 
groups, for which the Dirac method breaks down. Yet our method works, the modi- 
fied inner product being given by an expression similar to that in the compact group 
example above, with due care taken about the choice of the dense subspace L. 

To see how the covariant method operates in a system without constraints, con- 
sider an unparametrized relativistic particle, with Lagrangian L = —m\J\ — x 2 , 

where x (X . Ob . Ob 3 ) and x = dx./dt. The covariant phase space S cov of this model 

(Kijowski 1973) is 1R 8 , with co-ordinates {t,p, x l ,p 1 }. The covariant Legendre trans- 
form (Gotay 1991, Gotay et.al. 1993) leads to the primary constraint $ = p—H = 0, 
where H = ^p 2 + m 2 = u p . (Here H happens to coincide with the usual Hamil- 
tonian, but in field theory it would be a covariant generalization thereof.) The 
symplectic form on S cov is uo = —dp Adt + dp' 1 A dx l . The solutions to the equations 
of motion are precisely tangent to the null direction of i*u at any given point. The 
reduced phase space S® ov of S cov with respect to the constraint p — H is the usual 
phase space T*M 3 . The co-ordinates p and t have been eliminated by the reduction 
procedure. Time evolution on S° may be described by any observable acting as a 
Hamiltonian on S® ov , such as H. 

To quantize, we replace S cov by the covariant Hilbert space 7i cov = L 2 (M 4 ) (with 
■ip a function of (t,x)). According to our method, we have 



We find H° = L 2 (R 3 ; d 3 p/(2n) 3 ), and V defined by / p (t,x) = exp(-icu p t + ip ■ x). 
These functions are a complete set of generalized solutions of the constraints, which 
explains in part why alternative methods (e.g., Hajicek 1994) are equally well capable 
of handling this example. This is not generic, however (cf. the last example below). 
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4 The Polyakov particle 



To mimick minisuperspace, we now treat the relativistic particle in the parametrized 
Polyakov formalism (cf. Gotay et.al. 1993 for the classical part). The Lagrangian is 

where x = dx/dr, and g = goo is a metric on the real line. The covariant phase space 
is parametrized by {r,p, g, n}, subject to the primary constraints p — H = 0, 

with 

and 7r = 0. The symplectic form on S cov is u = —dp A dt — dp^ A dx^ + dn A dg. 
According to the general method of Gotay et.al. (1993) to handle parametrized 
theories, the secondary constraints are picked up by the covariant momentum map, 
in this case with respect to the action of Diff + (R), under which the Lagrangian is 
covariant. If N(r)d/dr is an arbitary element of the Lie algebra of this group (that 
is, a vector field on IR), this momentum map is given by & N = —pN — 2ngN, which 
is linear in the canonical momenta, as it should, and generates diffeomorphisms in 
the appropriate manner. 

The complete set of constraints, then, is p — H = ir = §n = (the latter for all 
N). In full general relativity, there would be analogous constraints, with g and 7r 
replaced by the lapse and shift components of the metric, and the diffeomorphism 
constraints now being of the form / s <7*<3>7v = (here a is a section of the covariant 
primary constraint bundle (Gotay et.al. 1993), and E is a Cauchy surface). In any 
case, the diffeomorphism constraints are linear in the canonical momenta also in 
that case. The primary constraints, i.e., ir = and p = H, are imposed on the 
finite- dimensional covariant phase space. It is only when p = H is substituted into 
the diffeomorphism constraints (and evaluated at a fixed time) that one quarter of 
the latter become quadratic in the canonical momenta, ruining their key property 
of generating four-dimensional diffeomorphisms, as well as the property that their 
Poisson algebra represents the Lie algebra of the group. Such a substitution would, 
in fact, lead to the ADM formalism. 
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In the covariant Hamiltonian formalism, on the other hand, the constraints are of 
the type we can handle with our quantization method, since the diffeomorphism con- 
straints are essentially the (covariant) momentum map, and the finite-dimensional 
constraint p = H, which is typical to the covariant method, and independent of 
the diffeomorphism invariance, can be thought of as a momentum map for an ac- 
tion of R. Also, the formalism clearly distinguishes between the lapse g and the 
infinitesimal diffeomorphism N{r)d/dr. 

In the present case, instead of reducing with respect to Diff + (IR), the same 
reduced phase space is obtained if we reduce with respect to an arbitrary one- 
parameter subgroup generated by some N, as long as iV > 0. Also, since g > 0, we 
can rewrite the constraints as $1 = p = 0, $ 2 = n = 0, and $ 3 = —^{p^—m 2 ) = 0. 
Having reduced by $i and $2, we then simply have to reduce T*R 4 by $3, which 
of course leads to the correct result for the reduced phase space; the restriction to 
p > has to be imposed by hand (cf. Landsman 1995). 

The quantization of S cov is H cov = L 2 (R 6 ), with ip a function of (r, x^, log(g)). 
The modified inner product is given by 




defined on ip, (p e L = C C (1R 6 ). Imposing p > 0, it follows that the physical Hilbert 
space H° and the map V are as in the previous treatment (the f p are now regarded 
as functions of (r, x M , log(g)), which happen to be independent of r and g). Again, 
one easily verifies that observables such as the generators of the Poincare group act 
in the correct way. 

5 A minisuperspace example 

We finally turn to a simple minisuperspace example (cf. Halliwell 1990). The La- 
grangian is 

£ = -^ 3 VV-0V^), 

which describes an FRW universe with radius exp(a), filled with a homogeneous 
massless scalar field (here k — 0, — 1 (open), or +1 (closed)). The covariant phase 
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space S cov is parametrized by {r, p, a,p a , 0,p</,, g, n}. The constraints are as in the 
previous example, except that the Hamiltonian is now given by 

H = ^e- 3a (-p 2 a +pl- Ke Aa ). 

After elimination of r and g, the physical phase space is obtained by reducing T*M 2 
by the constraint H = 0. 

In a suitable parametrization (of which the reduced space is evidently indepen- 
dent), which corresponds to putting $ 3 = ±(— p 2 a +p\ — nexp(Aa)), for k — 1 the flow 
generated by this constraint is given by 0(A) = 0(0) + p^X, P<f,(X) = p</,(0); a(X) = 
ilog[ v / 2^/cosh(2 v / 2E(A - A ))], p a (X) = v / 2^tanh(2 v / 2E(A - A )), where E = 
^(p a (0) 2 + K-exp(4a(0))) and A is determined by a(0) = \ \o^[\f2Ej cosh(2 v / 2EA )] 
(or p a (X ) = 0). Observables are, for example, the functions F 1 = p^ and 
F 2 = — (arctanh(p a /v / 2~E'))/2v / 2 _ E', which in fact project to canonical co- 
ordinates /i,/2 on the reduced space S°. These are observables of DeWitt type, 
e.g., F 2 (z) is the value of 0(^(A)) at the A for which p a (z(X)) = 0, z E T*R 2 . For 
k — —1 the above expression for the flow is formally correct as well, but in that case 
the classical motion is incomplete at infinity and a special interpretation is necessary 
(cf. Feinberg and Peleg 1995 for similar cases). 

The idea of taking the matter field (which is not an observable) as a time 
variable (cf. Isham 1993) is implemented in a parametrization-invariant way by using 
fi as the physical Hamiltonian of the model. If a (z) is the value of exp(a(^(A))) 
at the A for which 0(z(A)) = 0, then a physical prediction of the model would be 
the time evolution al(f 1 ,f 2 ,t) = fi/ cosh(2(/ 2 + t)). This describes a recollapsing 
universe, cf. Marolf 1995a. 

6 The wavefunction of the universe 

We quantize on H cov = L 2 (M 4 ), with -0 a function of (r, a, 0, log(g)). The mod- 
ified inner product is defined as in the previous example. For k — 0, 1, $3 is the 
unique self-adjoint extension of the usual Schrodinger quantization of $3 on the core 
C C °°(M 4 ). For k = the analysis is trivial, leading to H° = ©±L 2 (M; dk/2n), and 
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\/(°) = V}°> © V m defined through /^(t, a, 0, log(g)) = exp(ik(a ± 0)). For k = 1 
the physical Hilbert space is 7^5 = L 2 (JBL; dk/2ir), and some functional analysis (for 
which cf. Leis 1979, Picard 1989) shows that one has 

ff\r, a, 0,log(<?)) = 7r- 1 e^ V / s]nli(^^^| fe | /2 (Ie 2Q ). 

For k = — 1 the operator $3 is not essentially self-adjoint on C^°(1R 4 ); its deficiency 
indices are (1, 1), and one obtains a one-parameter family of self-adjoint extensions 
by specifying boundary conditions at infinity. For all of these, H°_ 1 = H®, and 
the wavefunction of the universe (see below) turns out not to be affected by the 
particular choice of the extension. For simplicity we choose the boundary condition 
specified in Picard (1989). This leads to 

if 1} (r, «, <f>, logfa)) = ie^cosech(7r|fc|/2)( J m/2 + J- m/2 ){\e 2a ). 

(Different boundary conditions at infinity would have led to a relative phase between 
Ji\k\/2 and J-i\k\/2-) 

These are solutions to the Wheeler-DeWitt equations of the model (which do 
not lie in Tt cov ), but the difficulty lies in the fact that now there are other linearly 
independent solutions as well. Our technique gives a procedure of choosing between 
them. In alternative quantization methods, the boundary conditions (at zero radius) 
singling out the relevant solution are lost. 

One may compare these solutions with proposed 'wavef unctions of the universe' 
\l> (cf. Halliwell 1990). This comparison is feasible if one has a rationale for letting 
k — > in fjf\ which might be provided if the physical Hamiltonian is chosen to 
be p^, as in the classical case. Up to a constant normalization factor, the Hartle- 
Hawking proposal here gives Io rather than our Kq for k — 1, and agrees with our 
Jo for k — — 1 (see Zhuk 1992). Vilenkin's wavefunction, which is always complex, 
differs from ours as well as Hartle-Hawking's in both cases. But note that according 
to Halliwell (1990) both proposals are somewhat ambiguous and not always well- 
defined - in fact, one can find contradictory statements in the literature concerning 
these proposals. Our expansion functions f^, on the other hand, only depend 
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on the domain chosen for the quantum constraints, which at least for k — 0, 1 is 
uniquely given (within the bounds of reason). 
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